Measuring protein concentration with entangled photons 
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Optical interferometry is amongst the most sensitive techniques for precision measurement. By 
increasing the light intensity a more precise measurement can usually be made. However, in some 
applications the sample is light sensitive. By using entangled states of light the same precision can 
be achieved with less exposure of the sample. This concept has been demonstrated in measurements 
of fixed, known optical components. Here we use two-photon entangled states to measure the 
concentration of the blood protein bovine serum albumin (BSA) in an aqueous buffer solution. We 
use an opto-fluidic device that couples a waveguide interferometer with a microfluidic channel. These 
results point the way to practical applications of quantum metrology to light sensitive samples. 



Even the most advanced sensors are bound by a hard 
limit in precision — the shot noise or standard quantum 
limit (SQL) that arises from statistical fluctuations. In a 
conventional optical interferometer for example the pre- 
cision with which an unknown optical phase <p can be 
measured is limited to 8<j) = 1/ \/N, where N is the (aver- 
age) number of photons used to probe (j) [T]-[6] . Increasing 
N is usually possible, by increasing laser power for ex- 
ample. However, in some scenarios the practical limits 
of laser power are reached and increasing the integra- 
tion time will reduce the bandwidth of the measurement 
below that required — gravity wave interferometers are a 
key example [7] . In other scenarios the sample to be mea- 
sured may be sensitive to light, such that one would like 
to minimise the photon flux or total number of photons 
that the sample is exposed to in order to reach the re- 
quired precision; put another way, one wishes to gain the 
maximum information allowable by the laws of physics 
for the given perturbation of the sample. It is this lat- 
ter scenario that we are focussed on here. By harnessing 
quantum superposition and entanglement the SQL can 
be overcome — quantum metrology enables the more fun- 
damental Heisenberg limit of precision, 5<j) = 1/N, to be 
reached pp. However, practical applications of quantum 
metrology require that samples of interest are integrated 
with quantum optical circuits. 

We use the optofluidic device shown in Fig. 1, con- 
sisting of a microfluidic channel that passes through one 
arm of a Mach-Zehnder interferometer (MZI), fabricated 
by femtosecond laser microfabrication [8HT0]. This de- 
vice combines the stability of integrated optics for high 
visibility quantum and classical interference [TTJ [12] with 
high precision handling of fluid samples [8] [13] . When a 
solution is fed into the microfluidic channel, any relative 
phase shift of light (and thereby concentration-dependent 
refractive index) in the sensing arm with respect to that 
acquired in the reference arm can be estimated from the 



interference fringes. 

The period of the interference fringes and the measure- 
ment precision 5<p depend on the particular state of light 
used to probe the sample. For classical states of light 
Scj) > 1/VN, the SQL. Going beyond this limit requires 
quantum states of light; the canonical example is the 
NOON state (|N0) + |0N))/>/2, which can achieve super- 
sensitivity and saturate the Heisenberg limit 5<j) = 1/N. 
Here, super-resolution results in a fringe periodicity that 
is 1/N times shorter than the one obtained with clas- 
sical light pp. The generation and detection of NOON 
states with large N is an active area of research [T4HT6] . 
We test the operation of our device for the N = 2 NOON 
state, which enables super-sensitity and super-resolution, 
and can be generated from two single photons input into 
a beamsplitter [17]. Photon losses in the interferome- 
ter reduce the measurement sensitivity, making NOON 
states non-optimal in general [18]. However, for two pho- 
tons, the NOON state is optimal [19] [20], an d we calcu- 
late that for the MZI configuration the best obtainable 
performance is achieved when input and output beam 
splitters have 50:50 splitting ratios (see Supplementary 
Information) . 

The device (Fig.[l]) was fabricated by femtosecond laser 
micromachining in a fused-silica sample to enable the 
integration of optical waveguides [9j [21] and microflu- 
idic channels [22] [23] in a three-dimensional architec- 
ture [10 : waveguides with slightly elliptical cross section 
were fabricated by astigmatic shaping of the writing laser 
beam [24 so that a small birefringence is induced to pre- 
serve linear polarization; the microchannel was fabricated 
by irradiating a double pyramidal structure followed by 
etching in a hydrofluoric acid solution in order to have 
perfectly vertical walls [22] [23] (see Supplementary Infor- 
mation). 

Photon pairs at A = 785 nm were generated via sponta- 
neous parametric down- conversion (SPDC) in a nonlinear 




FIG. 1. Quantum metrology in an optofluidic device, (a) Schematic of the experimental setup: A pump laser at X p — 392.5 nm 
generates pairs of downconverted photons at \ s — \% — 785 nm in a Bismuth Borate (BiBO) crystal. IF: interference filter, 
cl: collection lenses, PMF: polarization maintaining fibers, FA: fiber array, (b) Schematic of the Mach-Zehnder interferometer 
(MZI) interfaced to the microchannel. The fluidic channel has rectangular cross-section 500 /im x 55 /xm and extends from 
the top to the bottom surface of the glass substrate (~1 mm thickness). The MZI consists of two 50:50 directional couplers 
and has two arms of equal geometrical length; one waveguide crosses perpendicular to the microchannel, while the other 
passes externally, (c) Top image of the optical-fluidic interface, (d) Picture of the device with several interferometers and 
microchannels on chip, together with the fiber arrays for coupling input and output light. 



bismuth borate BiBaOe (BiBO) crystal (see Supplemen- 
tary Information) and collected into polarization main- 
taining fibres (Fig.JTJa)) . The photon pairs were coupled 
into the MZI via fiber arrays. Hong-Ou-Mandel (HOM) 
interference [25j [26] at the first directional coupler gen- 
erates a two photon NOON state [17 (|20) + |02))/>/2 
which is the state we use to probe the sample. Before 
interfering in the second directional coupler, the sensing 
mode acquires a relative phase shift in the microchannel 
that crosses the sensing arm of the MZI. 

Output photons are collected by an array of standard 
telecommunication fibres (monomodal at 1550 nm wave- 
length, but multimodal at the 785 nm wavelength used), 
in order to increase collection efficiency, and detected by 
four single-photon avalanche photodiodes (A, B, C and 
D in Fig.^a)), after non-deterministic separation at two 
50:50 fiber-splitters. With this detection scheme we are 
able to monitor the different two photon components of 
the output state |11), 1 20) and 1 02) and renormalize the 
measured fringes with respect to drifts during the mea- 
surement in the coupling between fiber arrays and MZI 
and source pair production rate (see Supplementary In- 
formation). 

To establish the quality of quantum interference in our 
device, we performed a HOM experiment [25], filling the 
microchannel with distilled water. A lossless MZI com- 
posed of two 50:50 couplers, with a relative phase shift 
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FIG. 2. Quantum interference in the Mach-Zehnder interfer- 
ometer when the microfluidic channel is filled with distilled 
water. The coincidences at the detectors A and B are plotted 
as a function of the relative delay between the two photons. 

<j) between the two arms, is equivalent to a beamsplit- 
ter of reflectivity cos 2 (0/2). However, asymmetric losses 
between the interferometer arms limits the maximum vis- 
ibility to: 

4T-(T-1) 2 

VHOM = ( T + 1 )2 W 

where T is the transmissivity across the microchannel, 
equivalent to the ratio of the transmissivity of the two 
arms and (j) = tt/2. 



3 




Concentration [%] 



FIG. 3. Quantum interference fringes. Normalized single pho- 
ton counts (a) and two photon coincidences (b) for different 
concentrations of bovine serum albumin in a buffer solution 
(full circles). The solid line represent a fitting of the exper- 
imental points with a sinusoidal curve. Error bars on data 
points are the same size of the dots and computed assuming 
Poissonian statistics of the detection events; other error, aris- 
ing for example from evaporation, are not taken into account. 

Figure [2] shows a typical two-photon detection rate 
across the two output modes as a function of the rel- 
ative arrival time, controlled with a translation stage 
(Fig. [IJa)). The visibility V = 86.7 ± 1.3% is almost 
ideal for a ratio of transmissivity T=61% since an upper 
bound of 88% is calculated from Eq. [I] We note that it 
is merely a coincidence that insertion of distilled water 
at room temperature in the channel results in <fi = tt/2 
to within the precision of this measurement. 

To test the operation of our device with a real sam- 
ple we chose Bovine Serum Albumin (BSA) in aqueous 
buffer solutions as a model fluidic sample that is stable 
and well-characterized [27]. Insertion of the solution in 
the microchannel was achieved by casting a droplet on 
the top aperture and exploiting spontaneous filling by 
capillary action. This geometry is chosen for its sim- 
plicity and extension to more sophisticated microfluidic 
delivery of the solution is straightforward [T3] . 

We performed sensing measurements using one photon 
and two photon inputs for 15 different concentrations of 
BSA ranging from 0% to 7% in 0.5% steps. Cleaning of 
the microchannel with deionized water and acetone was 
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FIG. 4. Refractive index change in the buffer solution as a 
function of BSA concentration. Experimental data (dots) are 
shown together with a linear fitting (solid curve). 

performed before and after each measurement. The sin- 
gle photon fringe is obtained by coupling only one pho- 
ton from the SPDC pair into the MZI and counting the 
number of detections from one output of the MZI. Figure 
|3ja) shows the single photon count rate normalized with 
respect to the sum of the singles from the two outputs 
(see Supplementary Information) together with theoret- 
ical fit function A/| 10 ) = (1 + Vi p / l cos(c/) + c/>o))/2, where 
Vi p h is the fringe visibility, (j) = aC (with C concen- 
tration of BSA and a constant), and c/>o is a constant 
phase offset term. A visibility Vi p / l =94±2.2% is esti- 
mated from the fit, compared to a theoretical prediction 

Vtheory 

,ip/i=97%, calculated taking into account the de- 
vice losses. 

Two photon fringes were measured by coupling photon 
pairs into the two input waveguides of the MZI and de- 
tecting coincidences from the two separate output chan- 
nels. Coincidence events C|n\ are normalized with re- 
spect to the sum of all the possible two photon outputs 
(see Supplementary Information) C|n)+C|2o>+C|o2> an d 
shown in Fig. [3^b) together with the theoretical fit func- 
tion 

AT (l + ^cos(20 + 0o)) 

7V |n> = 2 ' ' ' 

where the period is half that of the single photon 
fringe due to super-resolution. The visibility of the 
fit is V r 2 P / l =82±4.8%, in agreement with the theoret- 
ical prediction for the interferometer including losses 
Vtheory,2ph = 88%- This value exceeds the threshold for 
supersensitivity [28l|29]: V^ L = 70.7%. Current source 
and detector efficiencies prevent beating the standard 
quantum limit [30]. In addition to the 2.1dB (T = 61%) 
loss across the microfluidic channel we estimate 0.5 dB 
propagation losses and 0.5 dB bending losses, which fur- 
ther increase the visibility threshold. 
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The refractive index change An s of the BSA solution 
can be related to the phase shift 0, acquired by light 
during propagation in the sample, according to: 



An., 



A 



(3) 



where L = 55 \im is the microchannel length and 
A = 0.785 /im is the wavelength. The dependence of An s 
on the BSA concentration C (Fig.|4| can thus be inferred 
from the two-photon coincidences reported in Fig. [3J The 
experimental points are well fitted by a linear function, 
whose slope dn s /dC = 1.79 ±0.04 x 10 -3 is in very good 
agreement with the value of 1.82 x 10 -3 , previously re- 
ported [31 at A = 0.578 fim. 

Measurement of the concentration of a protein in solu- 
tion with entangled states in an integrated quantum pho- 
tonics device shows the potential for quantum interfero- 
metric measurement of light sensitive samples. Heralded 
or deterministic generation of larger entangled states will 
enable greater sensitivity [16 , when combined with high 
efficiency photon sources and detectors. Quantum optical 
circuits of that herald the generation of up to four pho- 
ton entangled states for quantum metrology have been 
demonstrated with lithographic waveguides [T5] . Multi- 
pass schemes [5] would be compatible with the optoflu- 
idic architecture demonstrated here, provided low loss 
switches can be integrated. More sophisticated microflu- 
idic delivery systems could be integrated for particular 
applications [8j [13]. Adding more waveguide capabili- 
ties, for example polarization-based quantum measure- 
ments [3j [32] , would enable measurement of samples that 
induce a concentration dependent rotation of the probing 
light polarization. 
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APPENDIX 

Fabrication of the optofluidic device: The inte- 
grated optofluidic device is fabricated using a Yb:KYW 
cavity- dumped mode-locked oscillator, which delivers 300 
fs pulses at 1030 nm wavelength, with 1 MHz repetition 
rate. The structures are fabricated by direct laser in- 
scription with sample translation enabled by high preci- 
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sion three-axes air bearing stages (Aerotech Fiber-Glide 
3D). Laser irradiation for both interferometers and mi- 
crochannels is performed in the same run. Several de- 
vices, each composed of a microchannel interfaced with 
a Mach-Zehnder interferometer, are realized on the same 
2 cm x 2 cm glass chip. 

Waveguides are fabricated using five overlapped scans 
at 1.5 mm/s with 300 nJ pulses at the second harmonic of 
the laser (515 nm wavelength), focused at 400 fim depth 
into the substrate by a 0.6 NA microscope objective. A 
cylindrical telescope, shaping the writing beam, is em- 
ployed to control the shape of the waveguide cross-section 

The microchannel is realized using the technique de- 
scribed by Vishnubhatla et al [23] . In order to obtain per- 
fect verticality of the microchannel walls we need to com- 
pensate for the intrinsic tapering induced by the etch- 
ing process along the z-axis (Figj5|. This is achieved by 
irradiating a helicoidal path in order to obtain double- 
pyramidal structure. A single laser scan at 0.5 mm/s, 
using 300 nJ pulses of the second harmonic, is employed, 
with the same focusing set-up used for the waveguides. 
After laser irradiation the sample is immersed in HF so- 
lution for 70 min. The HF acid selectively etches the 
irradiated material, proceeding along the z-axis from the 
bottom and top surfaces, simultaneously. This technique 
produces a microchannel through the 1-mm thick glass 
sample with a uniform and rectangular cross-section of 
55 fim x 500 fim. Polishing of the edges of the chip is 
finally performed, to allow effective fibre butt-coupling 
to the interferometers. 

Photon pair source: The reported experiments used 
degenerate 785 nm photon pairs, generated using spon- 
taneous parametric downconversion in a 2 mm thick, 
type-I nonlinear bismuth borate BiE>306 (BiBO) crys- 
tal phase matched for non-collinear generation (3° open- 
ing angle). The photon pairs were filtered using a 3nm 
full-width at half-maximum interference filter (IF) and 
focused with aspheric collecting lenses into polarisation 
maintaining fibres (see Fig. 1. A of the main text). The 
down conversion crystal was pumped by a ^ 10 mw 
392.5 nm beam (focused to a waist of ujq ~ 40 /im), gen- 
erated via collinear up-conversion — using a 2 mm thick 
BiBO crystal — of an attenuated 785 nm, 180 fs pulsed 
Ti:sapphire beam. The 392.5 nm light was filtered using 
four successive dichroic mirrors. 

Normalization procedure: Both classical and quan- 
tum optical interferometric experiments, such as the ones 
we performed on BSA solutions, generally consist in mea- 
suring a specific output state rate. Considering the 2 
x 2-ports Mach-Zehnder interferometer employed in our 
experiment (Fig. [lja)), classical measurements can be 
performed feeding the device with single photons from 
one input and counting output states 1 01 ) or 1 10) . Quan- 
tum measurements with two-photon NOON states can 




FIG. 5. Schematic of the femtosecond laser irradiated path for 
the fabrication of the microchannel. The double-pyramidal 
shape is adopted to obtain perfectly vertical walls after the 
etching process, which proceeds from the top and bottom sur- 
faces of the chip as indicated by the blue arrows. The final 
shape of the microchannel is depicted by the dashed box. 

be performed, on the other hand, feeding the interfer- 
ometer with two indistinguishable photons and detecting 
the output states |11), 1 20) or |02). The detection rate of 
any of these states depends on the interferometric phase 
(which is what we want to measure), as well as on pho- 
ton losses and detector efficiencies. If a single output 
state is monitored, the measurement result is prone to 
external perturbations, and variations in the losses or in 
the efficiencies cannot be distinguished from variations in 
the interferometric phase. Therefore we choose to adopt 
robust normalization procedures both for classical and 
two-photon quantum measurements, which correspond to 
single photons and coincidences detection, respectively. 

To build the single photon interference fringe, the fol- 
lowing counts are considered: 

^l|meas = P\ 10)^^1 

-42| meas = p\oi)T] A N2 

#l|meas = P|01)^B^1 
#2| meas = P\ 10 )T] B N2 

where Xk\ meas stays for experimental single count rates 
at detector X when photons are injected into input k; 
p\ij) is the probability to find a \ij) state at the inter- 
ferometer output, Nk is the rate at which input photons 
are effectively injected in the device and r]x is the effi- 
ciency with which a state \ij) is actually detected, and 
includes losses, fiber coupler splitting ratio and detector 
efficiency. The dependency from the unknown parame- 
ters, Nk and r]x, can be cancelled by retrieving p| i> from 
the experimental count rates according to the following 
formula: 

\/-^2 1 m eas B 1 1 meas / A \ 

Pl01} ~ / 401 ml 7 fJTT R9 | W 

V I meas- 01 1 meas ~T" \/ ^i-L | me as-° z |meas 
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FIG. 6. Schematic of the Mach-Zehnder interferometer. Lin- 
ear losses in the interferometer arms are modeled as extra 
beamsplitter (in red). 



The value of p| i> can thus be adopted as the normalized 
quantity for constructing the interference fringe. 

An analogous normalization method is employed for 
the evaluation of the two-photon interference fringe. Co- 
incidence counts at detectors AB, CD, AC and BD are 
acquired, which can be written as: 

AB\ meas = N^^A^B 
C£>| m eas = N^TJcVD 
^C|meas = N^VaVC ' 2 
BD\ meas = N^VBVD ' 2 



where XYl 



is the experimental coincidence count 



rate for the couple of detectors XY and is the rate 
of the \ij) state at the output of the interferometer. The 
factor 2 in the expressions of AC\ mediS and BD\ meas takes 
into account that the two photons can come from two 
equivalent paths. Also in this case we introduce a nor- 
malized quantity pin), which is independent from losses 
and efficiencies, denned as: 



P\ii) 



in) 



^lii) 



tot 



|n> 



■ N\ 



|20) 



■Ni 



|02) 



(5) 



where iV| 2 o> = ^02) > when a MZI with two 50:50 couplers 
is used. The former expression can be retrieved from the 
acquired coincidence counts according to: 



P|ii> 



sCDI 



(6) 



As detailed above, all the experimental quantities 
Xk\ meas or XF| meas depend on losses and efficiency 
terms, which in principle could be characterized in each 
set-up, but in real applications could anyway vary, e.g. 
due to misalignments induced by external perturbations, 
thus introducing unpredictable fluctuations in the mea- 
surements. On the contrary, the normalized quantities of 
Eq. [4] and |6] completely avoid this issue. 

Mach-Zehnder optimal configuration: For a Mach- 
Zehnder interferometer (MZI) with losses, the probability 




FIG. 7. Phase sensitivity S calculated from Eq. [15] as a 
function of the reflection coefficient of the two couplers. 



of a 2-photon coincidence can be expressed as [29] 

P llW = ,(i±^£i^±M, (7) 



where both the visibility V^-ph an d the coefficient rj de- 
pend on the losses in the two arms of the interferometer. 
From this equation we can compute the variance in the 
coincidences probability as 



AP n (<t>) 2 = P n (<j>)(l - P u (<t>)) 



(8) 



while the actual variance of the phase estimation is given 
by: 



l2 _ AP n O) s 



(9) 



In order to estimate the probability Pn (0) we analyse the 
MZI with losses in the Heisenberg picture. Linear losses 
in the interferometer arms are modelled as directional 
couplers where vacuum is injected (Figj6l) in one of the 
inputs. These ancillary modes are labelled and a\ 
while the real inputs of the MZI are a\ in and in . 

We will now calculate the probability of a coincidence 
event as a function of the reflectivity of the two di- 
rectional couplers of the interferometer, considering the 
losses in the two arms as parameters. Assuming that 
the first directional coupler has a reflectivity ri, with 



it will transform the creation operators of 



the input modes according to 



l,in 



jtia[ 



2, in 



(10) 



The supplementary couplers have transmissivity coeffi- 
cients that are related to the channel losses by Li = t\ . , 
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therefore the linear losses transform the creation opera- 
tors into 



r Ll a'l+jt Ll al iin 



(ii) 



The phase shifting operator acting on the sensing arm 



gives: 



d'"\ 

2 



d"\e il 
-J 



(12) 



while the second beam splitter generates the final output 
given by: 



a l,out 
a 2,out 



r 2 a">\+jt 2 a>"\ 
jt 2 a">\ + r 2 a'"\ 
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From Eq. [13] it is possible to calculate the probability of 



I 

a coincidence event as a function of the splitting ratio of 
the MZI directional couplers, which is given by: 



Pn(0) = |2r 1 t 1 r 2 t 2 (r| i e 2 ^ + r^)+r Ll r L2 (l-4r^ + 2(^+^))e- 



(14) 



Once the coincidence probability is known it is possible 
to calculate the maximum sensitivity S of the device for 
a two photon measurement given by the square root of 
the inverse of the phase variance: 



^ dPn(0) ^2 

S 2 = max^dcj) 2 )- 1 = max <f>( A p^^ 2 )' ( 15 ) 



I 

Figure [7] shows the phase sensitivity S as function of the 
reflection coefficients Ri = r\ of the two directional cou- 
plers of the interferometer with a maximum for the case 
where R\ = R 2 = 0.5. 



